Equivalence of interest rate models and lattice gases 



(N 

o 

(N 

Oh- 
< 



PL, 
U 

d 

<4-H ' 
I ' 



> 

in 

OS 

o 
o 



Dan Pirjol 1 

J J. P. Morgan, 277 Park Avenue, New York, NY 10172 

We consider the class of short rate interest rate models for which the short rate is proportional to 
the exponential of a Gaussian Markov process x(t) in the terminal measure r(t) — a(t)exp(x(t)). 
These models include the Black, Derman, Toy and Black, Karasinski models in the terminal 
measure. We show that such interest rate models are equivalent with lattice gases with attrac- 
tive two-body interaction V(ti,ta) = — Cov (x(ti), xfo)). We consider in some detail the Black, 
Karasinski model with x(t) an Ornstein, Uhlenbeck process, and show that it is similar with a 
lattice gas model considered by Kac and Helfand, with attractive long-range two-body interactions 
V(x,y) = — a(e i\ x ~ y \ — e An explicit solution for the model is given as a sum over the 

states of the lattice gas, which is used to show that the model has a phase transition similar to that 
found previously in the Black, Derman, Toy model in the terminal measure. 
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I. INTRODUCTION 

We consider in this paper the class of one-factor inter- 
est rate models with log-normally distributed short rate 
in the terminal measure. In these models the short rate is 
driven by one Gaussian Markov process x(t). Such a pro- 
cess is defined by two conditions: i) for any set of times 
t\ < t-2 < ■ ■ ■ < tk, the values x(t 2 ), ■ ■ ■ ,x(tk)) 

have a joint normal distribution; ii) the evolution of x(s) 
for all s > t depends only on x{t). It can be shown that 
the most general process of this type is a time-changed 
Brownian motion, and includes the Ornstein-Uhlenbeck 
process as a particular case 

This class of models includes the Black, Derman, To; 
(BDT) model, and the Black-Karasinski (BK) 
model, formulated in the terminal measure. The terminal 
measure is sometimes used in practice for these models 
[H, as opposed to the spot measure in which the models 
were originally formulated, due to the ease of calibration 
and simulation. Such models have been also proposed as 
approximations to the Libor market model [5|, [6| , and as 
particular parametric realizations of Markov functional 
models 0, 0] ■ A choice of measure amounts to a distri- 
butional assumption for the dynamical variables of the 
model. See [8| for a readable introduction to the related 
concepts of martingales and measure for stochastic pro- 
cesses, and their relation to arbitrage pricing theory. 

In this paper we show that these interest rate models 
are equivalent with lattice gases with attractive two-body 
interaction V(ti,t2) — —Cov(x(ti),x(t2)), placed in an 
external potential. The solution of the models can be 
expressed explicitly as an expression for the one-step zero 
coupon bond given by a sum over occupation numbers in 
the lattice gas. The expectation values required for the 
simulation of the model correspond to thermodynamical 
potentials in the lattice gas model. 

We discuss in some detail the Black, Karasinski model 
with constant mean reversion 7, which is equivalent to a 



lattice gas with attractive two-body interaction V(x, y) = 
— a(e _7 ' 2:_! '' — e~ llyX+y ^). This is similar to a lattice gas 
model considered by Kac H, Kac, Uhlenbeck, Hemmer 
[l(]| and Kac, Helfand (lUll^. This model generalizes 
the BDT model in the terminal measure, which corre- 
sponds to 7 = 0, and is equivalent with a Coulomb lat- 
tice gas with attractive two-body interactions. The latter 
model was studied in Ref. [llj], where it was shown that 
it displays discontinuous behaviour in volatility, which is 
similar to a phase transition in condensed matter physics 

[HE!. 

The equivalence with the lattice gas models suggests 
alternative simulation methods for these interest rate 
models, which express expectation values as sums over 
the states of the lattice gas. For small lattices this can 
be done by explicit summation over the lattice gas states, 
while for bigger lattices efficient numerical methods are 
available from statistical mechanics, such as Gibbs sam- 
pling and the Metropolis algorithm. We illustrate this 
approach by a numerical study of the BK model, which 
shows that the volatility phase transition observed in the 
BDT model in Ref. [l3[ persists also for this model. 



II. THE INTEREST RATE MODEL 

We consider a short rate interest rate model in discrete 
time. The model is defined on a finite set of dates 



= t < U < 



< U, 



(1) 



For simplicity we will assume that U are equally spaced, 
and denote r = ti+i — t% with i = 0, 1, • • • , n — 1. 

The fundamental dynamical quantities of the model 
are the zero coupon bonds Pij = P^u ■ They are defined 
as the price at time t% of a payment of 1 made at time 
tj. They are stochastic quantities, and can be expressed 
as functions of an one-dimensional Markov process x(t). 
For definiteness we consider in the following that x(t) is 



2 



an Ornstein-Uhlenbeck process with zero mean reversion 
level 



dx{t) = -~fx(t)dt + adWlt) . 



(2) 



The mean and variance of x(t) conditional on x(0) = 
are 



[x(t)\x(0) = 0] = 



(3) 



E[x 2 (t)\x(0) = 0] = _(1 - e-^) = G(t) . (4) 

The arguments of this paper can be easily extended to the 
more general case of x(t) an arbitrary Gaussian Markov 
process. By the Doob's representation, the most general 
Gaussian Markov process can be represented as a time- 
modified Brownian motion [l[ 



(5) 



*(*) = /(*) / 9(s)dW(s) 



with f(t),g(t) deterministic functions of time, and W(t) 
a Brownian motion. 

We define the Libor rate (or simply Libor) for the 
(fcj, ij+i) period as 



Li = T~ 



1 



1 



(6) 



The model is defined by specifying the functional depen- 
dence of the Libor rate Li on the Markov driver x(U) 



Li exp 



(x{u) - hs{u) 



(7) 



where Li are constants to be chosen such that the initial 
yield curve Pqj is correctly reproduced. This implies that 
the Libors Li are log-normally distributed in the terminal 
measure. 

This model is similar with the Black-Karasinski model 
@, S] , up to the difference that the latter is usually for- 
mulated in the risk-neutral measure, while in the model 
considered here the short rate Li is expressed in terms of 
x{t) defined in the terminal measure. 

In the limit when the time step is taken to zero r — > 0, 
this model becomes a continuous time short rate model, 
and the short rate r(t) — lim T _>o L t / T (t) satisfies the 
stochastic differential equation 



dr(t) 
r(t) 



(a(t) - j In r(t))dt + <rdW{t) 



(8) 



with a(t) a function depending on Li and a. We recognize 
this as the short rate evolution in the Black-Karasinski 
model Q. 



is a martingale. The mathematical statement of this re- 
sult is expressed as 



V(t)/N(t) =E[V(T)/N(T)\F t ], 



(9) 



for any t < T. This holds under fairly general assump- 
tions, among which market completeness is the most im- 
portant one. Speaking loosely this means that the model 
contains sufficiently many tradeable instruments to allow 
any possible payout to be reproduced as a combination 
thereof. 

The choice of the numeraire N{t) is not unique, 
and any particular choice defines a measure for the 
stochastic process followed by the discounted asset prices 
V(t)/N(t). Two particular choices are most common in 
the context discussed here. The spot measure, or the 
risk-neutral measure, takes N(t) to be the money mar- 
ket account at time t, while the terminal measure (or 
i n -forward measure) takes N(t) — P ti „ to be the zero 
coupon bond maturing at time t n . Once the condition 
([7|1 is imposed, different measure choices produce differ- 
ent observable distributional properties of the dynamical 
quantities of the model (rates and bonds), and thus ef- 
fectively correspond to different models. 

We will work in the terminal measure in the following. 
It is convenient to introduce the zero coupon bond prices 
divided by the numeraire Pt, n , which will be denoted as 
Pi.j = Pi,j/Pi, n - They are martingales in the terminal 
measure, and thus satisfy the condition (J9j> , which reads 
explicitly 



Pi. 



E 



P, 



P 



k.n 



(10) 



for all i < k < j < n. The one-step discounted zero 
bond Pi^+i(xi) will play an important role in writing 
the analytical solution of this model. It satisfies a few 
conditions, following from the martingale condition (|f 01) . 
First, its expectation value is known in terms of the initial 
yield curve 



E[P iii+1 (xi)} =P ,i+i. 
It also satisfies the two conditions 



(11) 



Pi,i+i(xi) = E[P j;+ i, l+2 (xi + i)(l + L t+ i(x i+ i)T)\Fi] 

(12) 

P 0}i =E[P i>i+1 (x i )(l + L i (x i )r)] (13) 



The first condition (|T2|) determines recursively the func- 
tional form of Pi 7 i + i(xi), starting with P n ~i. n = 1 and 
proceeding backwards in time. This is given explicitly as 
a conditional expectation value 



A. Explicit solution of the model 

According to the fundamental theorem of arbitrage 
pricing theory Q, the price of a financial asset V(t) ex- 
pressed in units of a simpler asset N(t) (called numeraire) 



P, l+1 ( Xl ) = E[ J] (l + £*re**-* G *)|.F<]. (14) 

k=i+l 

The second condition (|T3"]) can be used to determine 
Li also recursively, once Pi.i+i{xi) has been determined, 



3 



using the relation 

~ Li = , P ^- P W . (15) 

E[P M+i exp(xi - ^Gi)]r 

For simplicity we denote the value of the Markov driver 
at time as cCj = x(ti), and its variance as G(ti) = G{. 



We will state in the following the closed form of the 
solution of this model. The solution expresses the dis- 
counted one-step zero coupon bonds P^i+i as a sum 
of terms containing 0, 1, 2, • • • , n— i — 1 Lj factors. Writ- 
ing the first few terms explicitly this is given by 



Pi, i+1 (xi) = 1 + £ L jT exp{w i ~ i Xi ~ 2 w2 °" l)G «) 
j=i+i 

n— l 

+ £ L J L k T 2 exp((w^+w k -*)x l --(w^+w k -*) 2 G l + X^ +■■■ (16) 

j>k=i+l 

+ £ £ L h L h ■ ■ ■ L Jk r k exp ( £ u^"**, - ±(£ ^"*) 2 G 4 + £ X jadb ) . 

k<n-i-lS k eTi a=l a=l l<a<b<k 

I 



We denoted here the weight w = exp(— jt), and the auto- 
covariance of the Markov process x(t) as 

X jk = Cov(x(tj),x(t k )) (17) 

2 

= f_( e -7l*j-tfcl _ e —r(tj+tk)\ 
2 7 

The general term in Eq. (116[) containing fc < n — £ — 1 
factors of Lj is given by a sum over all subsets S k = 
{jiiJ^r ■ ■ ,jk} of fe indices chosen from the n — i — 1 
indices 7* = {£ + 1, £ + 2, • • • , n — 1}. 

In the limit of zero mean reversion 7 — > 0, we have 
w = 1 and G(t) = cr 2 t, and the expression (ITB1) simplifies 
drastically. In this limit all terms with the same number 
of Lj factors have the same functional dependence of Xi , 
and we recover the simple form obtained in Ref. [l3| 

n-l 

A,i+i(i) = £c5V**-^ aG < (is) 

i=o 

where the coefficients are given by 

4 i) =£%4- 1 ---^ fc exp( £ A^-J (19) 

S k Ka<b<k 

where Xjjt = cr 2 min(ij,ifc). In [lj| these coefficients 
were determined recursively from a recursion relation, 
see Eq. (12) in Ref. [la ]. Equation (fl9| gives an explicit 
solution of this recursion relation. 

An important role is played in this model by the ex- 
pectation values of the form 

A^)=E[P M+1 e^-^ 2G *] (20) 

n-l 

= 1 + £ Lj-rexp^u^Gs) H 

j=i+i 



+ £ £ ^31^32 ■ ■ ■ Lj k r k 
k<n—i—l Sk 
k 

x exp^G* £ u^> ~* + £ X ja Jb ) . 

a=l l<a<b<k 

We enumerate in the following the applications of these 
expectation values with <f> = 0, 1, • • • . 

The expectation value of Pi.i+i (corresponding to <p = 
0) is constrained by the requirement that the initial yield 
curve Po,i is correctly reproduced, see (fTTjl . 

n-l 

E[P M+1 ] = fj),<+i = 1 + £ L jT + • • • (21) 

j=i+l 

-£/•.; /•; ■'' exp( £ + 

S fc Ka<b<k 

The sum on the right-hand side is linear in Li+\ and thus 
can be used to solve explicitly for this constant, provided 
that all Lj with j = i + 2, • • • , n — 1 are already known. 
This is given in Eq. (1151) in a form more convenient for 
practical calculation. 

The (f> — I expectation value appears in the calculation 
of the convexity-adjusted Libors Li Eq. (fT5|) . which can 
be written equivalently as 

Finally, A^(j) with j E > 1 determines the j— th 

moment of the Libor distribution in its natural (forward) 
measure according to the relation [ll| 

E«+i[(£<) j ] - ^— (ii^lPM+i^* 4- ^ '] 

-P 0,i+l 
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Pa 



-(L t y e -^-^N t (j) (23) 



In the limit of zero mean-reversion 7 — > the above 
expectation values are given by simple expressions [l3[ 



J G,1 



3=0 



(24) 



For sufficiently small volatility <r, the expectation val- 
ues Ni(<f>) given in Eq. (|20|) can be computed in an ex- 
pansion of the small parameter L{T <C 1, and keeping 
only the terms linear in this parameter is sufficient for 
most applications. In this approximation we have 



Ni(cf>) 



E 

j=i+i 



L^ d Te 4 



0((4 wd r) 2 ) 



(25) 



The distribution of the Libors in their natural measure is 
approximatively log-normal and the ATM caplet volatil- 
ity is 



2 

CT LN 



U 



(26) 



In the model with zero mean reversion 7 = 0, it was 
noted in Ref. [l3| that for volatility a above some crit- 
ical value, the higher order terms in the expansion (|20[) 
become comparable to the linear terms of O(Lj-r). The 
actual expansion parameter becomes L^r exp(cr 2 ^) and 
terms of all orders in L^ wd r become important. This 
leads to a discontinuity in the hrst derivative of the ex- 
pectation value Ni((f>) with respect to the volatility a, 
which is similar to a phase transition in condensed mat- 
ter physics fl4l[l5l|. 

In the next section we express the expectation values 
(f20| as averages over the grand canonical ensemble in an 
equivalent lattice gas model. This is used to show the 
existence of a phase transition also in this model, using 
a numerical simulation. 



B. Proof 

The result (fT6)) can be proven using the following basic 
identity. For any numbers n k = 0, 1 associated with the 
ordered sequence of times t = t < t\ < t 2 - ■ ■ < tjy, 
the following expectation value with x(t) the Ornstein- 
Uhlcnbeck process ((2]) is given by 



E 



exp(^n fc (a; fe - -Gfe))|.Fj 

fc=i 

N ^ N 

exp (x t J2n k e^ tk - -Gt^nke" 7 **) 

k—l 

1 N 



(27) 



fc=i 



lt k \2 



where Xj^ is the covariance of the process x(t) given 
above in Eq. (fT7|) . This is a slight generalization of an 
identity used in Ref. 0,[ll| to compute the partition func- 
tion of a lattice gas with exponential interaction. It can 
be easily generalized to the case of a general Gaussian 
Markov process x(t). 

The discounted one-step bond Pi^i+i (xi) is given by the 
conditional expectation (|14p. Expanding out the product 
yields terms with 0, 1, 2, • • • factors of L k T 1 up to n — i — 
1 factors. There are ( n ~^~ 1 ) terms containing TV such 
factors, and they are given by a sum over all subsets 
{nk} = {nkx , nk 2 , • • • , rik N } of N indices out of the total 
of n — i — 1 indices. A generic term has the form 

N 

£ nf =1 (L fci r)E[exp (5>*>*, - ^G^.))^] 

{n k } 3 = 1 

= ]Tllf =1 (L fcj T) (28) 
{»*} 

iV 1 N 

fe=l k=l 

+ E X k j ,k l n k j n kl 

kj<ki 

where the expectation value was computed using the 
identity (|27l) . This reproduces the terms containing N 
factors of Z^r in Eq. (fli?|) . This completes the proof of 
(HI. 



III. THE LATTICE GAS MODEL 

The interest rate model considered in the previous sec- 
tion is equivalent with a one-dimensional lattice gas with 
attractive long-range potential 



V{x,y) = -a{e-^ x - y \ - e - l{x+v) ) 



(29) 



The particles of the lattice gas are constrained to sit at 
positions X4 = ri, with i = 1, 2, • • • , n— 1. The n sites of 
the lattice gas are labeled as j = 0, 1, • • • , n— 1. The sites 
j are in one-to-one correspondence with the discrete set of 
simulation times {tj} of the interest rate model. At each 
site at most one particle can be present. We define nj the 
occupation number of the site j. It can take values or 
1, depending on whether the site j is vacant or occupied. 
The Hamiltonian of the lattice gas model is 



H = Y,< 



^ Sjknjnk 



j=l 3>k=l 

The two-body interaction is 

-7t|j— k\ _ „-~/T(j+k) 



£ : jk: 



a ( e -7T|J-*| _ e -7TU+KJ) 



j^k=l 



and the single-site energies are 



(30) 



(31) 



(32) 
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For the application to the interest rate model we are 
interested not only in the entire lattice system, but also 
in the subsystem Tj of the lattice consisting of the sites 
Tj : {i + 1, • • • , n — 1}, in total nj = n — i — 1 sites. 

Assume that the subsystem Tj of the lattice gas is 
placed in a position-dependent chemical potential 

H®(t)=iiG i e-v<t- t *'> (33) 

The grand partition function of the subsystem Tj of 
the lattice gas with the Hamiltonian (|3T)|) and placed in 
the chemical potential (|3"3"|) is given by 

n—i—X 

zdp,T)= E cx p(-' 9i7 +' 3 E ^ ( %))( 34 ) 

n=o s N jes N 

The sum over the number of particles A~ runs from to 
n — i — 1, the number of lattice sites in the subsystem 
Tj. For each A~ the sum runs over all configurations Sn 
of A" occupied sites, which are subsets of A^ sites of the 
n — i — 1 sites in the system Tj. 

The correspondence of this lattice gas model with the 
interest rate model is realized through the following re- 
lation between the grand partition function -2-j(/x, T) and 
the expectation value (|20|) 

Ni(4>) = Zi(n,T) , (35) 

provided that the parameters of the lattice gas are related 
to those of the interest rate model as 

2 

a/3 = °— , 4> = fo (36) 

This system is similar to the one-dimensional gas con- 
sidered by Kac and by Kac, Uhlenbeck, Hemmer 
A lattice version of the gas model, very similar to that 
considered here, was examined by Kac and Helfand in 
Ref. [ill EH • More precisely, the latter papers consider 
a lattice gas, where the particles occupy a lattice with 
A" nodes and lattice spacing 1, and interact by two-body 
attractive potentials V(\x — y\) — —wye *> x ~y\. This 
model has a phase transition in the so-called van der 
Waals limit, which is obtained by first going to the ther- 
modynamical limit of large N, followed by the infinite 
range limit 7 — > 0. In the van der Waals limit the lattice 
gas model has a liquid-gas phase transition with critical 
temperature (3 c a = b, and the equation of state is given 
by the van der Waals equation supplemented by the equal 
area rule 10]. 

At this point it may be useful to recall a few well-known 
facts about phase transitions in one-dimensional systems 
17]. Although a phase transition does not exist in a one- 
dimensional system with short range interactions [18( , it 
is possible for such a system to have a phase transition 
provided that the interaction is sufficiently long range. 
Sufficient conditions which have to be satisfied by the 
interaction in order for a phase transition to exist in a 
one-dimensional system were given in [191 ] . The papers 
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FIG. 1: (Color online) Plots of In Ni(l) vs a for several values 
of the mean-reversion parameter 7, with i = 30 in a simulation 
with n = 40 quarterly time steps r = 0.25. The black curve 
(leftmost) corresponds to 7 = and is obtained using the 
method used in Ref. [131] . The other curves (from left to 
right) are obtained by explicit summation over the occupation 
numbers of the lattice gas as explained in the text: 7 = 0.1% 
(blue), 1% (red), 2% (green), 5% (orange). 

[loj provided the first instance of phase transition in a 
one-dimensional system, and showed explicitly that this 
can occur in a system with long-range interactions. The 
results of [l(| have been extended to more general inter- 
actions and higher dimensional systems in Ref. (20j . 

The zero mean-reversion limit of the interest rate 
model 7 = is the Black, Derman, Toy model in the ter- 
minal measure [lj|, and is equivalent with a lattice gas 
model with attractive Coulomb two-body interactions, 
placed into an external potential. This can be seen by 
writing the covariance of the Markov driver for this case 
as 

-V(ti,t 2 ) = Cov(x(ti),x(t 2 )) y =o = cr 2 min(ii,i 2 ) 
= \v\\h-t 2 \-{h+t 2 )). (37) 

The first term describes an attractive linear interaction 
between the pair of particles at sites t\,t 2 , while the sec- 
ond term can be represented as their interactions with 
the repulsive external field of a static charge placed at 
the site i = 0. 

The one-dimensional gas with Coulomb interaction be- 
tween several types of charges was studied, using meth- 
ods very similar to those employed here, by Edwards and 
Lenard [2l|. Our Coulomb lattice gas is different from a 
usual Coulomb gas in that all particles attract each other. 
The thermodynamics of a one-dimensional system with 
linear attractive potentials was considered in Ref. , al- 
though periodic boundary conditions were imposed such 
that the resulting form of the interaction is different from 
that considered here. A connection between stochastic 
processes and the (two-dimensional) Coulomb gas was 
realized in a different context in Ref. 23] . 

The lattice gas with non-zero mean-reversion consid- 
ered here differs from that studied by Kac and Helfand 
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[Tl],[l2j in several respects, due to the peculiarities of the 
interest rate model. 

1. The presence of the Lj factors requires the introduc- 
tion of single-site energies £j associated with the lattice 
sites. These energies are different and thus the space ho- 
mogeneity of the system is lost. This space homogeneity 
was crucial for the analytical solution of the model in 
the thermodynamical limit . A similar approach is 
unlikely in this case for this reason. 

The single-site energies Sj are constrained by the condi- 
tion (|22|) such that the initial yield curve Po,i is correctly 
reproduced. According to this relation, Sj depends on 
the properties of the subsystem Tj-% of the lattice gas, 
and must be determined by a recursive procedure start- 
ing with the smallest subsystem T n _2 and adding one 
lattice site at a time. 

2. The two-body interaction in the lattice gas ([2T?]) con- 
tains a second exponential term exp(— jiU+tj)), which is 
not present in Rcfs. [HI [HI- This is due to the fact that 
the expectation values (|2"T)) are conditional on x(0) = 0, 
while [Til [l2| integrate over x(0). While the new term 
does not have the typical form of a two-body interaction, 
its inclusion does not present any problem of principle. 
Also, this term becomes vanishingly small if the subsys- 
tem Tj is chosen such that 3> 1, and the simple ex- 
ponential Kac interaction is recovered in this limit. 

The equivalence of these interest rate models with lat- 
tice gases suggests an alternative way of calibrating and 
simulating such models. The expectation values iVj(</>) 
are usually [J, 0] computed by evaluating the nested in- 
tegrations over the values of the Markov driver x(t) at 
the simulation times, using numerical approaches such 
as finite difference or Monte Carlo methods. The results 
([TO]) and ([2"0"]l suggest that the expectation values Ni(<j>) 
can be also computed as averages over the grand canon- 
ical ensemble in the lattice gas. For small lattices, this 
can be done by explicit summation over all possible oc- 
cupation numbers (2 n configurations for a lattice with n 
sites), while for larger lattices alternative methods famil- 
iar from statistical mechanics can be used, such as Gibbs 
sampling and the Metropolis- Hastings algorithm [U, [25[ . 

As an illustration of this approach, we show in Fig. [1] 
the results of a simulation of the BK model in the termi- 
nal measure performed by summing over the occupation 
numbers of the lattice gas. These plots show the mul- 
tiplicative convexity adjustment ln-/Vj(l) for i = 30 as 
function of a for several values of the mean-reversion pa- 
rameter 7. The simulation assumed n = 40 quarterly 
time steps r = 0.25, for a total simulation time t n = 10 
years. The forward yield curve is flat with Lf" = 5%. 
The 7 = curve is obtained using the recurrence method 
of [l3| . and the remaining curves were obtained by com- 
puting Ni(l) using (|20p by explicit summation over the 
2«-»-i _ states of the subsystem T30 of the lattice 
gas. 

These results show that the transition observed in 
Ref. [IH persists also in the model considered here. The 
mean-reversion 7 allows one to control the range of the 



two-body interaction in the lattice gas. In the 7 — !• 
limit the lattice model particles attract each other with 
Coulomb potentials, while for 7 7^ the potential be- 
comes exponential and is given in Eq. (j2U)) . In the 7 — > 
limit the results of [13J are recovered: the convexity 
adjustment factor increases suddenly above the critical 
volatility a CI ~ 32%. As the mean reversion 7 is in- 
creased from zero, the transition persists, and the critical 
volatility increases from its 7 = value. The 7 — > limit 
is well-behaved, as expected for a finite size lattice. 

The study of the 7 = limit of this model presented 
in Ref. [H[ showed that the phase transition is not visi- 
ble under usual simulation methods used in practice for 
such interest rates models, such as finite difference or 
Monte Carlo methods. This is due to the fact that these 
methods effectively truncate the range of values of the 
Markovian driver x(t) to a few (~ 5) multiples of a\ft. 
Such a truncation omits the contributions to the expec- 
tation values Ni(4>) which are responsible for the phase 
transition. The alternative method proposed here offers 
a possible way to study the properties of these models, 
free of these limitations. 



IV. CONCLUSIONS 

We presented in this paper the exact solution of a class 
of interest rates models with log-normally distributed 
short rates in the terminal measure. The solution is for- 
mulated naturally in terms of a lattice gas with sites cor- 
responding to the simulation times of the model t t . At 
each site only one particle can be present, and the parti- 
cles interact by attractive two-body potentials Vij which 
are determined by the stochastic process followed by the 
short rate. 

The analogy with the lattice gas models simplifies very 
much the simulation of these models, as many of the 
important expectation values in the interest rate model 
can be written in closed form as averages over the grand 
canonical ensemble in the corresponding lattice gas. The 
numerical evaluation of these averages is straightforward 
for small lattices (few simulation times in the interest 
rate model), while for larger lattices the number of con- 
figurations (2" for a lattice with n sites) becomes too 
large for direct evaluation, and approximation methods 
familiar from statistical mechanics may have to be used 

We used the exact lattice gas solution to study numer- 
ically the Black, Karasinski model in the terminal mea- 
sure with constant mean-reversion and volatility. This 
showed the appearance of a phase transition in the con- 
vexity adjustments of single-period interest rates, similar 
to that noted in the Black, Derman, Toy model in the 
terminal measure in Ref. [13| . This adds further support 
to the suggestion made in Ref. [l3[ that the presence of 
such a transition is generic for all interest rate models 
with log-normally distributed rates in the terminal mea- 
sure. Although the present numerical study considered 
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only the version of the model with constant parameters, 
the method can be extended without any major difficulty 
also to the more general case of time-dependent model 
parameters. This is in contrast to the method of the re- 
cursion relations used in to solve the 7 = limit 
of the model with uniform volatility, which does not ap- 
pear to be easily extended beyond this case due to the 
unmanageable complexity of the resulting expressions. 

The equivalence of the interest rates models consid- 
ered with interacting lattice gases shows that the former 
have a rich dynamics which has not been fully explored. 
Physical intuition about the lattice gas equivalent should 
give further insight into the dynamics of the interest rate 
models. In particular, one natural question is whether 



a phase transition similar to that studied in Ref. [10| is 
present also in the lattice gas considered here, and if it is 
observed also for a finite size lattice. The analog of the 
van der Waals limit for this case corresponds to simul- 
taneously scaling the volatility as a — 007 as the mean 
reversion is taken to zero 7—7-0. It would be interest- 
ing to see if the behaviour of the system in this limit 
has implications also for the practically relevant case of 
non-zero volatility. 

Finally, it would be interesting to investigate whether 
the exact solution presented here can be extended also 
to other interest rate models, with more general distribu- 
tional properties. Hopefully the lattice gas analogy will 
remain useful also for more general interest rate models. 



[1] J. L. Doob, Heuristic approach to the Kolmogorov- 

Smirnov theorems, Ann. Math. Stat. 20, 393 (1949). 
[2] F. Black, E. Derman and W. Toy, A one-factor model 

of interest rates and its application to Treasury bond 

options, Fin. Anal. J. 46(1), 33 (1990). 
[3] F. Black and W. Karasinski, Bond and option pricing 

when short rates are lognormal, Fin. Anal. J. 47(4), 52 

(1991). 

[4] L. Andersen and V. Piterbarg, Interest rate modeling, 
Atlantic Financial Press, 2010. 

[5] A. Daniluk and D. Gatarek, A fully lognormal Libor mar- 
ket model, Risk, 2005 

[6] O. Kurbanmuradov, K. Sabelfeld and J. Schoenmak- 
ers, Lognormal approximations to Libor market models, 
J. Comp. Finance 6(1), 69, 2002. 

[7] J. B. Hunt and J. E. Kennedy, Financial derivatives 
in theory and practice, Wiley Series in Probability and 
Statistics, 2005. 

[8] M. Baxter and A. Rennie, Financial calculus: An in- 
troduction to derivative pricing, Cambridge University 
Press, 1996. 

[9] M. Kac, On the partition function of a one-dimensional 

gas, Phys. Fluids 2, 8 (1959). 
[10] M. Kac, G. E. Uhlenbeck, P. C. Hemmer, On the van 

der Waals Theory of the Vapor-Liquid Equilibrium. I. 

Discussion of a One-Dimensional Model, J. Math. Phys. 

4, 216 (1963); ibid. 4, 229; ibid. 5, 60. 
[11] M. Kac and E. Helfand, Study of several lattice systems 

with long-range forces, J. Math. Phys. 4, 1078 (1963). 
[12] E. Helfand, Approach to a phase transition in a one- 
dimensional system, J. Math. Phys. 5, 127 (1964). 
[13] D. Pirjol, Phase transition in a log-normal Markov 

functional mode l, J. Math. Phys. 52, 013301 (2011), 

arXiv:1007.0691 



[14] E. H. Stanley, Introduction to phase transitions and crit- 
ical phenomena, Oxford University Press, 1987. 

[15] N. Goldenfeld, Lectures on phase transitions and the 
renormalization group, Frontiers in Physics, Addison- 
Wesley, 1992. 

[16] D. Pirjol, Nonanalytical behaviour in a log-normal 
Markov functional model, larXiv: 1104.0322 1 

[17] D. C. Mattis, The many-body problem: an encyclopedia 
of exactly solved models in one dimension, World Scien- 
tific, New York 1992. 

[18] L. Landau and E. M. Lifschitz, Statistical physics 1, Perg- 
amon, Oxford, 1980. 

[19] F. Dyson, Existence of a phase transition in a one- 
dimensional Ising ferromagnet, Comm. Math. Phys. 12, 
91 (1969); ibid. 21, 269 (1971). 

[20] J. L. Lebowitz and O. Penrose, Rigorous Treatment of 
the van der Waals-Maxwell Theory of the Liquid- Vapor 
Transition, J. Math. Phys. 7, 98 (1966). 

[21] S. F. Edwards and A. Lenard, Exact statistical mechanics 
of a one-dimensional system with Coulomb forces. II. The 
method of functional integration, J. Math. Phys. 3, 778 
(1962). 

[22] A. Isihara, Consideration of a phase transition in a one- 
dimensional gas, Physica 64, 497 (1973). 

[23] F. Dyson, A Brownian motion model for the eigenvalues 
of a random matrix, J. Math. Phys. 3, 1191 (1962). 

[24] N. Metropolis et ai, Equations of state calculations 
by fast computing machines, J. Chem. Phys. 21, 1087 
(1953). 

[25] W. K. Hastings, Monte Carlo sampling methods using 
Markov chains and their applications, Biometrika 57, 97 
(1970). 



